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ing semiclassical methods. Our approach is based on the Landauer formalism and the Fisher-Lee relations, 
appropriately generalized to spin-orbit interaction, and a semiclassical representation of Green functions. We 
calculate conductance coefficients by exploiting ergodicity and mixing of suitably combined classical spin-orbit 
dynamics, and making use of the Sieber-Richter method and its most recent extensions. That way we obtain 
weak anti-localization and confirm previous results obtained in the symplectic ensemble of Random Matrix 
Theory. 
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I. INTRODUCTION 

Ballistic transport through chaotic cavities realized as quan- 
tum dots in semiconductor heterostructures has been a cen- 
tral issue in mesoscopic physics for many years. The uni- 
versal transport properties observed in this context can be de- 
scribed on a phenomenological level by random matrix the- 
ory' (RMT). The same applies to disordered systems, where 
averages over impurities can be shown to be equivalent to 
random matrix averages. This not being possible for indi- 
vidual, clean cavities, theoretical explanations of the RMT- 
connection have been provided making use of semiclassical 
methods, which are based on the Landauer formalism^ and 
semiclassical representations of Green functions. This ap- 
proach^ leads to questions that are closely analogous to prob- 
lems arising in semiclassical explanations of universal spec- 
tral correlations in classically chaotic quantum systems. Re- 
cent progress in the latter context is based on the seminal work 
of Sieber and Richter* and its extensions^'^-^. This method has 
been adapted**'^''" to be able to successfully explain conduc- 
tance coefficients, including the effect of weak localization, 
i.e., a decrease of conductance at zero magnetic field. Fur- 
ther studies have been devoted to analyses of the universality 
of conductance fluctuations and of shot nois o'^'^^i'^ . 
(For an overview see, e.g., Ref. 10). 

In the work mentioned transport properties were consid- 
ered for ballistic, non-relativistic electrons, neglecting their 
spin. In the emerging field of semiconductor based spin elec- 
tronics'^ (spintronics), however, one requires an efficient con- 
trol of the spin dynamics associated with electrons in non- 
magnetic semiconductors. This purpose calls for an inclu- 
sion of spin-orbit interactions into studies of transport proper- 
ties. In contrast to previous theories neglecting the spin, here 
one would expect appropriate classical spin-orbit dynamics to 
produce weak anti-localization, i.e., an enhancement of the 
conductance at zero magnetic field. This prediction is also 
obtained on the phenomenological level provided by RMT, 
where a half-integer spin requires the symplectic, as opposed 
to the orthogonal, circular ensemble. On this ground one ex- 
pects universal conductance fluctuations and shot noise also 



to be affected by the presence of spin-orbit interaction&i^i^. 
A first semiclassical approach'^ to these questions employs 
the semiclassical representation of the Green function in spin- 
orbit coupling systems derived in Ref. [Isl and considers the 
first order of the semiclassical Sieber-Richter expansion. It, 
moreover, assumes a randomization of spin states, which is 
shown to be responsible for weak anti-localization. 

In this paper our goal is to extend the results of Ref. [nitoall 
orders of the Sieber-Richter expansion, and to base the semi- 
classical estimates entirely on dynamical properties of suit- 
ably combined classical spin-orbit dynamics'^. These then 
replace the randomization hypothesis of spin states made in 
the analytic part of Ref. [l3- In order to determine the spin 
contribution to transmission amplitudes we closely follow an 
analogous calculation introduced in the context of semiclas- 
sical explanations of spectral correlations in quantum graphs 
with spin-orbit couplings^"-^' . We also comment on shot noise 
and on the variance of conductance fluctuations. 

As our model we consider a two dimensional cavity with 
two straight, semi-infinite leads with hard walls. Apart from 
boundary reflections, particles with mass m, charge e, and spin 
s move freely within the leads and are subjected to a magnetic 
field and to spin-orbit interactions inside the cavity. Although 
the relevant case of electrons enforces the spin to be = 1 /2, 
we deliberately allow for general spin s. Below this will allow 
us to point out characteristic differences between integer and 
half-integer spin. The Hamiltonian governing the dynamics in 
the cavity reads 

^=:^(p--A(x))'+s-C(x,p) . (1) 

Im \ c / 

Here A is the vector potential for an external magnetic field 
and C contains all couplings of the translational degrees of 
freedom to the spin operator s. This may include Zeeman- 
, spin-orbit, Rashba-, or Dresselhaus-type couplings. More- 
over, in order to model the hard walls we require Dirichlet 
conditions at the boundaries of the cavity and of the leads. 

The paper is organized as follows: Section II is devoted to 
a generalization of the Landauer formalism and the Fisher- 
Lee relations to systems with spin-orbit interaction. Then we 
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present semiclassical representations of S-matrix elements in 
that case. In Section III we first introduce ergodicity and mix- 
ing conditions that include a classical spin-orbit interaction. 
This is followed by our calculation of the conductance in two 
ways: in the configuration-space and in the phase-space ap- 
proach. In Sections IV and V we then outline how our ap- 
proach can be extended to calculate shot noise and conduc- 
tance fluctuations, respectively. An Appendix contains a cal- 
culation whose result is central to the phase-space approach 
employed in Section III. 



n. PRELIMINARIES 

We follow the usual approach to obtain semiclassical ap- 
proximations to transmission by employing the Landauer 
formalism- and introducing semiclassical representations for 
Green functions. In the absence of spin-orbit interactions this 
procedure is well established^^'^^'^"*. Here we briefly describe 
the extensions required by the presence of spin-orbit interac- 
tions (see also Ref.flTl). 



A. Landauer formalism with spin 

The Landauer formalism provides a link between conduc- 
tance coefficients, as defined through 



4 = D 



(2) 



and S'-matrix elements. In (|2]i the indices label the leads, V,,, is 

the voltage applied at lead m and /„ is the current through lead 

n. Here the number of leads may be arbitrary. An S'-matrix 

element 5"™ / is defined as the transition amplitude between 

an asymptotic incoming state in the lead m, characterized by 

the collection a,'„ of its quantum numbers, to an asymptotic 

outgoing state in the lead n, accordingly characterized by a„ . 
I II — I 

In Refs. 23'24 the Landauer formalism was derived from 
the Schrodinger equation in linear response theory, making 
use of an appropriate Kubo-Greenwood formula. We first re- 
mark that an inclusion of spin, interacting with the transla- 
tional degrees of freedom via a Zeeman, spin-orbit, Rashba, 
or Dresselhaus coupling, into this method causes no problems. 
Although the current density is modified, its conservation in 
the form required for the Kubo-Greenwood expression of the 
conductivity to hold is indeed guaranteed. On then obtains for 
transmission (i.e. m ^ «), 



and for reflection (i.e. m = n). 



gnn - ^ rdEf^{E) ( {2s + \)K - I 



(3) 



(4) 



Here N„ is the number of open channels in the lead « (without 
spin degeneracy) at energy E, and fp{E) denotes the Fermi 



distribution function at inverse temperature p. Of course, this 
requires the spin quantum number s to be half-integer. 

In a next step S'-matrix elements have to be related to Green 
functions G{x,x' ,E). These satisfy the equations 



(27;i(p-fA(x))^ 



-s-C(x,p) 

:5(x-x') 



e)g{x,x',E) 



and 



(5L(p' + £A(x'))'-£)g(x,x',£) 
+C*(x',p')G(x,x',£)s = 8(x-x') . 



(5) 



(6) 



The unusual form of the second equation is dictated by the 
fact that G(x, x' , £■) is a hermitian {2s + 1 ) x {2s + 1 ) matrix in 
spin space. In the following we will always choose advanced 
Green functions, fully characterized by Eqs. (|5]l and (|6]l as 
well as the condition that, asymptotically in the leads, they 
contain only outgoing contributions. 

As in the case without spin^^'^"* one can then express the 
5-matrix elements in terms of the (advanced) Green function. 
Up to a global phase factor, for m ^ n this yields 



and for m — n 




W„ rW„, 

djn I dy'^ sin 
h) 

Gao' (Xn , X„j , £ ) , 



(7) 



2h^ ^KM, f^-' 



X sm 



W„ JO 

w„ 



djn I dy'„ sin 



Goa'(x«,x^,£')+5o 



(8) 



Here we have introduced coordinates x„ = {x„^y„), where 
x„ > is a longitudinal, outward running coordinate in the 
lead n and < y,, < W„ is the corresponding transversal co- 
ordinate (see also Figure [T]l. The transversal quantum num- 
ber is Qn = 1,...,N„ with associated wave number ka = 



2mE /h" — a^K^ /W^. The number A^„ of open transversal 
channels then is the largest integer a„ that leaves the wave 
number real. Moreover, o = ~s, is a spin index such that 
altogether a„ = {E,a„,a). 

We remark that in Eqs. (|7J and ^ the points x„,xj„ can 
be chosen anywhere in the respective leads. For later con- 
venience we take them on the connection of the leads to the 
cavity, i.e., with x„ = = xj„. 



B. Semiclassical Green function and transmission amplitudes 

In order to proceed further, one requires a semiclassical rep- 
resentation for the Green function defined in Eqs. (|5]l and 
In Ref. 18 this was achieved through an asymptotic expansion 
in powers of Planck's constant h for the quantum propaga- 
tor generated by the Hamiltonian ([T]) which yielded, after a 
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I 



Figure 1 : Sketch of the geometry 



Fourier transformation, a respective semiclassical expansion 
for the Green function. The range of vaUdity of this proce- 
dure follows from the observation that, since the spin oper- 
ator s is linear in h, the energy scale of the spin-orbit inter- 
action term becomes small as compared to the kinetic term 
in the limit h ^ 0. This condition is equivalent to the spin- 
precession length being large compared to the Fermi wave- 
length. In semiconductor heterostructures this requirement is 
usually fulfilled. 

The semiclassical representation for the Green function ob- 
tained in Ref . [H reads 

G(x,x',£)~ Ay{x,x\E)exp{ii/h)Syix,x',E)) , (9) 

y(x.x') 

as /; — > 0. The sum extends over all classical trajectories 
y(x,x') generated by the classical Hamiltonian 



1 / e \2 
^^o(x,p) = — (p--A(x) 

2m \ c 



(10) 



(plus reflections from hard walls) that run from x' to x at en- 
ergy E. Choosing (x,x') — (x„,xj„) as in Q, the rele- 
vant trajectories are those that enter the cavity at lead m and 
leave through lead n. Moreover, Sy{x,x\E) is the classical 
action of the trajectory, and the leading order of the amphtude 
Ay(x,x'.E) reads 



Ay{x,x',E) = 



ih\/2%iTi 



fCyDy{x:,9',t){\+0{h)) . (11) 



Here Vy is a Maslov index of the trajectory y, and 









Cy . - — 


det ^ 


dxdE 1 









(12) 



The contribution of the spin is, in leading semiclassical order, 
completely contained in the spin-transport matrix DY(x',p',f). 
This is the spin-i representation of the spin propagator 
dy{x' jp^t), which is defined as a solution of the equation 

^dy{x',p',t) + '-C (X(0,P(0) •afl'Y(x',p',f) = (13) 



with initial condition c/y(x',p',0) = 1. Here (X(f),P(f )) is the 
point in phase space of the classical trajectory y at time t. Its 
initial point at time f = is (x',p'). Moreover, <J is the vector 
of Pauli spin matrices. Therefore, dy is an SU(2)-matrix that 
can be seen as a propagator for the spin along the classical 
trajectory y. 

Upon dividing the trajectory y into two pieces yi and y2, 
such that f = f 1 + f2, the spin propagator is clearly multiplica- 
tive. Since Dy arises from a group representation it inherits 
this multiplicative property from the propagator, i.e., 

DY(x',p',fi+f2)=DY2(X(ri),P(fi),f2)£'Y,(x',p',fi) . (14) 

This relation will be used extensively in Section Hill 

In order to obtain a semiclassical representation of trans- 
mission amplitudes we insert the expression (|9]l into equation 
(l7]l. Then the integrals over y and / , respectively, are evalu- 
ated, asymptotically as /i — > 0, with the method of stationary 
phase. In this context we stress the following important obser- 
vation: The number of accessible transversal states (including 
spin) in the n-th lead is (25+l)A^„ = (2i+ l)[\/2m£W„/(7iri)], 
where [x] denotes the integer part of x € K. We choose the 
widths Wn of the leads to formally shrink proportionally to Ti 
in this limit (compare also Ref. 14) and hence set W,, = W„/i, 
to the effect that the sin-factors in Eqs. (|7]l and (O contribute 
rapidly oscillating phases. These have to be taken into account 
when determining stationary points of the total phases in the 
integrals. The condition of stationary phase hence imposes the 
following restrictions on the transversal momenta. 



Py = 



and 



Py- 



dSy a„K 



(15) 



(16) 



upon entry and exit, respectively, of the trajectories. If the 
points of entry and exit are free of magnetic fields, and thus 
p = mx at these points, one can characterize the trajectories 
in terms of the angles 9 and 9', under which they enter and 
leave the cavity with respect to the longitudinal directions of 
the leads (see also Figure |2]). These angles are related to the 
transversal momenta ( fTSI l and ( fT6] l through sin 9 = py/ylmE 
and sin 9' = p\.l\f7mE. If one wished to keep the widths of 
the openings fixed, however, the method of stationary phase 
would enforce the conditions p\. = = py upon the trajecto- 
ries, thus leading to different semiclassical expressions than 
the ones we use henceforth. 

Collecting now all terms that emerge in the stationary phase 
calculation finally leads to the following leading semiclassical 
contribution to the S'-matrix elements. 



Ca;„~ L %e.e')'DY(e,e/) exp ((///!) 5^(9,6': 
Y(e,e') 



(17) 



where the sum extends over all trajectories that run from lead 
m through the cavity to lead n and are characterized by the 
conditions (flSl l, ( fTSI l, expressed in terms of the angles of entry 
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and exit. The explicit form of the factor ^^(e e') is the same as 
if there were no spin preseni22, 



III. SEMICLASSICAL CALCULATION OF 
CONDUCTIVITY COEFFICIENTS 



B. 



inh sgn(±a^,)sgn(±a„) 



2W,„W„ |cosecose'M2ig9,^|i/2 



X exp m 



(18) 



Here M^^g g,j is an element of the monodromy matrix of 
Y(9,9') that arises from the matrix appearing in ( fTSl l by a re- 
striction to the phase space directions transversal to the tra- 
jectory. Furthermore, /^^(e.e') is a modified Maslov index that 
contains the index Vy(e_e') from Eq. (fTTl i and additional phases 
resulting from the stationary phase calculation of the integrals 
overy,, andyj,,. 

The above result ( fTTb primarily refers to transmission am- 
plitudes (n 7^ m), but can be carried over to the case of reflec- 
tion (« = m). The reason for this is that the additional term 
^a„a'„ in ® is canceled by the contribution of direct trajecto- 
ries in the opening of the lead n that never enter the cavity^. 

The ultimate goal being a semiclassical calculation of the 
conductance coefficients ^ and (|4|i, one therefore requires 
the evaluation of double sums 



Y(9,e')y(9.e') (19) 
xexp((//;i) (Sy-Sy)) 



over classical trajectories. This will be the task for the rest of 
this paper 

To simplify the calculations, from now on we restrict our 
attention to the case of two leads. With an incoming wave in 
the lead m = 1 we are thus dealing with the transmission co- 
efficient g2i and the reflection coefficient gn. To this end we 
will determine the transmission matrix 5^' and the reflection 
matrix 5'', leading to the transmission and reflection coeffi- 
cient 



T 



n21 

a2a'. 



?22 



(20) 



oi2.ce. 



a2,Ol2 



respectively. Hence, at zero temperature the current through 
lead 2 is 



h = -{lVy + {^-{2s+\)N2)V2) , (21) 

h 



where T and are taken at the Fermi energy Ep- Together 
with the condition g2\ + g22 = 0, expressing that equal volt- 
ages at both leads produce no cuiTent, this yields the relation 



l2^-T (Vi-Vz) 
n 



(22) 



The calculation of the double sum iT% over classical tra- 
jectories requires input from dynamical properties of the asso- 
ciated classical system. With spin-orbit interactions present, 
one therefore first has to identify an appropriate classical sys- 
tem. Moreover, ergodic properties of the classical system im- 
ply necessary ingredients for the further calculation. The di- 
agonal contribution to the double sum is evaluated with a sum 
rule^ '^, whereas the non-diagonal terms are evaluated follow- 
ing the Sieber-Richter metho d'^'^i'^'^^ . 



A. Classical spin-orbit dynamics 

The classical dynamics that enter the semiclassical repre- 
sentation (|9]l consist of two parts^^: the motion of the point 
particle generated by the Hamiltonian ( fTOl i. including elas- 
tic reflections from hard walls, and the spin that is driven 
by this motion according to (fTsl l. These contributions can 
be combined into a single dynamics on a spin-orbit phase 
space^'^. The relevant classical trajectory is (X(f),P(f),g(f )), 
with initial condition (x',p',g') at f = 0. Here g G SU(2) and 
g{t) — dy{x' ,t)g provides the spin part of the combined mo- 
tion. We remark that this description of spin appears quantum 
mechanical. However, by passing to expectation values of the 
spin operator dy ^Ody in normalized spin states % (Heisenberg 

picture), the spin variable becomes a unit vector {%, dy jOdy%) . 
Hence the spin part of the combined phase space is a unit 
sphere. The two views of the spin motion, either on SU(2) 
or on a unit sphere, are in fact equivalents^. In both cases we 
will therefore speak of classical spin-orbit dynamics. 

Ergodicity is a concept developed for closed systems. It 
can, however, be suitably extended to open systems of the kind 
under consideration here. To this end one divides the config- 
uration space Q of the device into a closed part Qc, consisting 
of the cavity with the leads truncated and the openings closed, 
plus the infinite leads. From now on we suppose the shape of 
the closed part to form a chaotic billiard, ensuring ergodicity 
of the motion inside the cavity. Then p(f) is the probability 
for a typical trajectory to stay within the cavity at least up to 
time t. For large times. 



P(0 ^exp(-f/x) , f ^oo , 

with inverse dwell time 

I h , 

- = —r{Ni+N2) , 
T inA 



(23) 



(24) 



in which A denotes the area of the closed part Q^. For the 
associated part of phase space we also introduce the volume 

d^x [ d^p5{E-Hoix,p))^2nmA (25) 

of the energy shell. This expression has no integration over 
the spin part, since the Hamiltonian is independent thereof. 
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and an integration over SU(2) with respect to Haar measure 
dg yields one. 

For the open system the concept of ergodicity has to be 
modifed in that the possibility of a trajectory to leave the 
cavity must be taken into account. When the motion inside 
the cavity is ergodic this leads to the following relation be- 
tween phase-space averages and time averages over typical 
spin-orbit trajectories. 



lQj^xJ^2d^pJsu(2)dgfi'>i,P,g)5{E-Ho{x,p)) , (26) 



as r oo. Here / is an arbitrary function on the combined 
phase space, and (...) denotes an average over initial condi- 
tions. This relation, which properly reflects the chaotic nature 
of the combined classical spin-orbit motion, provides the basis 
for the further use of dynamical properties in the calculation 
of the sum ( fT9b over classical trajectories. 



The stronger mixing property, which we also assume to 
hold henceforth, means that correlations of two observables 
/ and h decay, i.e.. 



lim,^oe/g^y2^/R2 d'^pS^^j^^^dgh{X{t),V{t),g{t))f{x,Tp,g)h{E-HQ{x,T^)) 

r 



id^P!su{2) dg /!(x,p,g)5(£ -//o(x,p)) /g^ <iV/jj2 d^p' Ssu{2) dg'f{'ii',p',g')?>iE-Ho(x',p')) 



(27) 



B. Transmission and reflection coefflcients in the 
configuration-space approach 

In a first step we calculate the leading semiclassical contri- 
bution to transmission and reflection coefficients from equa- 
tion ( fTgl l, averaged over a small energy window, by using the 
configuration-space approach. Such a calculation has been 
performed previouslyii, however, with a sum rule that only 
takes the particle motion into account. The spin contribution 
was built in subsequently, assuming that traces of products of 
spin-transport matrices can be replaced by certain averages. 
Here we reproduce the result obtained in Ref. Tt* by using a 
sum rule for the complete spin-orbit dynamics that follows 
from ( l26l ). Thus we base the assumptions made in Ref. ITtI on 
a firm dynamical ground. 

As h ^ the terms in the double sum ( fT9] l are highly os- 
cillatory, except for contributions with 5y = 5y. Generically, 
if no symmetries are present, this only occurs for the diagonal 
y = y. In the event that time-reversal invariance is not bro- 
ken, however, the time-reversed trajectory y^^ has the same 
action as y. Of course, y^' is only among the trajectories to 
be summed over in the case of reflection (n = 1 = m) when, 
moreover, = 9'; i.e., only for 5^' / with ai — a', . All further 

terms are oscillatory, with a decreasing importance of their 
contribution, after averaging over an energy window, when 
the action differences increase. Below we calculate the two 
leading contributions to the quantity 



%i 

(28) 

for systems with time-reversal invariance: (i) the diagonal 
contribution in which the sum over Y is restricted to Y = y 
(for transmission) or y = y*' (for reflection), and (ii) the one- 
loop contribution in which the sums over y and Vare confined 
to so-called Sieber-Richter pairs (see also Ref.UTl). 



Due to the unitarity of the spin-transport matrices, in the 
diagonal case terms with Y = y yield a spin contribution of 
lLx{DyDlj) = 2i+ 1. Thus, the diagonal contribution to (|28] | 
can immediately be obtained from the respective result with- 
out spii 



25+1 



diag/A£ Ni + N2 



(29) 



fl J an additional 



In the case of reflection (n—l—m) with a i 
diagonal contribution arises from the terms with Y = y^\ if 
time-reversal invariance is unbroken. Its spin contribution is 
Tr(DyD^[) — Tr(Dy). One hence requires a suitable sum 
rule that incorporates the combined classical spin-orbit mo- 
tion. For this purpose we choose the function 

f{X{t),P{t),g{t)) 

= -5m 9) 8(x(r)) {@{y{t)) &{y{t) - W,)) (30) 
m 

xTr{KMt)g{0)-')f 

in (|26] |. Here Ks{g) denotes the spin-i representation of g G 
SU(2), "& is the angular variable in planar polar coordinates 
for p and &{y) is a Heavyside step function. An evaluation 
of ( |26] | with the function ( l30b then leads to the sum rule (as 



L |B/Tr(D2) 

y.Ty<T 



2Wi 2%mA Jo 



dtp{t) 



(31) 



After an average over a small window in energy this, together 
with ( |29] l, finally yields the semiclassical result 



E 

0,0' = -s 



aid 



diag/ AE 



2.+ l + (-l)2-'5„,„; 
N1+N2 



(32) 
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Figure 2: A Sieber-Richter pair of trajectories 



for the diagonal contribution to (l28T l. For 5=1/2 the right- 
hand side is l/(A^i +A^2)- 

Sieber-Richter pairs of trajectories are characterized by the 
fact that one trajectory possesses a self-crossing with a small 
crossing angle e, thus forming a loop. The partner trajectory 
then looks like the former one cut open at the self-crossing, 
but with the loop direction reversed and then glued together, 
such that the self-crossing is replaced by an almost-crossing, 
see Figure ID In principle, the trajectories in such pairs can 
have an arbitrary number of self-crossings, but the magni- 
tude of their contributions to (l28T l decreases with increasing 
numbers of places in which self-crossings are paired with 
almost-crossings. The most important ('one-loop') contribu- 
tion comes from pairs which differ in one crossing. In order 
to calculate the one-loop contribution one requires the distri- 
bution of the crossing angles e for pairs of trajectories with 
loops of duration T , 

Ps{£,T) = / d^x' [ d^p' r dtips{E,T,ti) . 

(33) 

Here ps(e, Tji) is a density of crossing angles defined as 

ps{£j,ti) = r " dts\J\5{E-HQ{P{ts))) 
Jo 

xTr(7r4g(f)(g(0))-'])'5(8-K(f.,f/)) 
x5(X(f,)-X(f, + f/)) , (34) 

where K{ts,ti) denotes the angle between the velocities v(fs) 
and v(fj +ti). Given a crossing angle e, the minimal duration 
for a loop to close is 7;,„„(e). In chaotic systems this quantity 
behaves like r„„„(e) = (9(loge) as e ^ 0^. Furthermore, 



|7| = |v(f,)xv(f, + f/)| 

|v(f,)| |v(f, + f/)| sinK(f,,f/) 



(35) 



is a Jacobian, and f,, ?/ denote the time along the trajectory up 
to the starting point of the loop and along the loop, respec- 
tively. 



Assuming that the classical spin-orbit dynamics are not 
only ergodic, but also mixing, the distribution ( |33T l can be cal- 
culated further. It can be identified as the left-hand side of an 
appropriate relation of the type (l27b . The right-hand side then 
yields, as e — > 0, 



{-ly 2E 
TzA m 



sme 



j-2 

TT . 

2 -* -'-miu 



(£)- 



T^- (e) 

mm V / 



(36) 

This expression differs from the respective one without spin 
that was obtained in Ref. 4 only by a factor (— 1)^'^ i.e., a 
sign in the case s — 1/2. With this information at hand the 
one-loop contribution can be calculated as in the case without 
spin^ 



, finally yielding 



.21 



o.o'=— J 



1— loop/ AE 



{Ni +N2) 



2 ■ 



(37) 



This is in accordance with what has been obtained in Ref. 



C. Transmission coefficients in the phase-space approach 

Higher orders in the 'loop-expansion' described above 
have been calculated previously for spectral form factors^ as 
well as for conductance coefficients for systems without spin 
contributions^. The approach taken in these papers utilizes 
trajectories in classical phase space and identifies the pairs of 
self-crossings/almost-crossings in configuration space as pairs 
of trajectories with almost-crossings in phase space, which 
differ in the way they are connected at the (almost) crossings. 
This point of view opens the possibility for a classification of 
the trajectory pairs in terms of their encounters^. Here we fol- 
low this phase-space approach and amend the previous result^ 
with the contribution of the spin-orbit coupling. 

To be more precise, we consider trajectories that possess 
close self-encounters (in phase space), in which two or more 
short stretches of the trajectory are almost identical, possibly 
up to time reversal. These stretches are connected by long 
parts of the trajectory, which we call loops. We then form 
pairs (Y,y) of such trajectories in which y and "/ are almost 
identical (up to time reversal) along the loops, but differ from 
each other in the way the loops are connected in the encounter 
region. In order to quantify these encounters, we introduce 
a vector v, whose l-th component, v/, denotes the number of 
encounters with I stretches. Hence the total number of en- 
counters is y = L/>2 vi, with a total of L = Lz>2 stretches 
involved. In general, however, given a vector v, there will be 
N{v) > I different trajectory pairs associated with it. These 
may, e.g., differ in the order the loops connect the encounters, 
or in the relative directions, in which the encounter-stretches 
are traversed. 

To reveal the phase-space structure of trajectory pairs and 
to compute their contributions to ( fT9l ) one introduces Poincare 
sections, which cut the trajectories into pieces. In order to 
adapt this cutting to the sequence of encounters and loops one 
chooses a Poincare section in every of the V given encounter 
regions. We then denote hy t'^j, j = l,...,la,Oi— 1 , • • • , V the 
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times at which the encounter stretches pierce this section, and 
by f^^., the duration of the encounters. To this cutting of the 
trajectories corresponds the splitting 



:Dl+iDl...Di 



(38) 



of the spin-transport matrices which, with an obvious nota- 
tion, follows from the composition rule (fT4l) . The spin trans- 
port along the partner trajectory then reads 



(39) 



Here r\j = ±1, depending on the relative orientation of the 
trajectory between the j — 1-st and the j-th cutting of y and 
respectively, through the Poincare section. We notice that at 
this point time-reversal invariance enters crucially. Moreover, 
the indices kj take care of the fact that in y and y' the loops may 
be traversed in different successions. Thus the spin-dependent 
weights in (1281 ) for each pair of trajectories are approximately 
given by 



Tr(Z)yZ)y) 



Tr(DL. 



.D 



tr|L^ 



(40) 



The calculation of transmission amplitudes performed in 
Ref. 1^ has now to be modified in that the expressions (|40] | 
must be included. To this end we recall the strategy devised 
in Refs. I^S For each encounter one introduces coordinates 
on the Poincare section adapted to the piercing by the trajec- 
tories and the linear stability of the dynamics. In encounter a 
the coordinates (i", m"), j = 1 1, describe the separa- 
tion of the j + 1-st piercing from the y'-th one along the stable 
and unstable manifolds, respectively, of the latter The total 
of L — y stable and unstable coordinates are then collected in 
the vectors (s,u). In these coordinates action differences of 
partner trajectories (approximately) read as 



A5 = Sy — 



(41) 



Moreover, the requirement that encounters be close can then 
be expressed in terms of the condition , |m"| < c with some 
constant c, which yields the duration of an encounter 



X max/ { I Si I } max j { | | } 

spin / 



(42) 



One then introduces a density (s,u) of encounters, 
weighted with the spin contribution, for trajectories of dura- 
tion T with a given encounter structure specified by the vector 
V. In analogy to the case without spin'^ this leads to the fol- 
lowing approximation. 



EE^(^) r ■■■ r d^'^ud^-^sexp{{i/h)AS) 



2\ 



(43) 



to the quantity 



nnd 



nil 
020', 



AE 



2s+l 



A£ Nl +N2 



(44) 



After summing over all possible values of 02, a[, this yields 
the non-diagonal contribution to the energy-averaged trans- 
mission amplitude T , compare ( l20t . ( |29] ). 

The essential point now is to calculate the density 
w^™(s,u). In the case without spin-orbit interaction the cor- 
responding expression wj- (s, u) was defined in Ref. 6 as a den- 
sity of phase-space separations s and u similar to the density 
P(e, r) with respect to e in the configuration-space approach. 
It was given as 



wr (s,u) = 



X rn«'o0(r-x:/„c-Eo) 

■'0 j=l \ a=l j=\ ) 



1 



/a 



xJl^ n8((X(^«.)>P(^a;))-^a.) 

a=l '<?nc \j=l 



(45) 



The average in the first line is over all possible initial points 
of the trajectory. In the second line the integration extends 
over all loop durations f/; their lengths are constrained by the 
theta function. In order to prevent over-counting^, the product 
of all encounter durations f^^. is divided out. The last prod- 
uct guarantees that the position of the orbit at times when it 
pierces through the sections are fixed as z^j- This denotes the 
first point of the orbit in which it pierces through a certain 
section plus the separation thereof as specified by the coordi- 
nates s and u. From Eq. ( l45T l one obtains WjP'"(s,u) by includ- 
ing Tr(DYDy) under the integral. Using that the durations of 
encounters are semiclassically large, compare (l42l l. the result 
can be obtained in analogy to ( |34] | by employing ( |27] |. The 
right-hand side then yields 



spin 



(s,u) 



(46) 



i.e. a factorization into the spin-independent part identical to 
W7(s,u) and a spin contribution 



M. 



SU{2) 

xTr 



5£/(2) 

TlTt 



(47) 



In order to calculate (l47T i we follow the method developed in 
Refs. I20II21I for the spectral form factor of quantum graphs 
with spin-orbit interaction. In analogy to Theorem 6.1 of 
Ref.|2l| we find in the present context that 



My^ = (2i+l) 



2i 



(-1) 

2s + 1 



L-V 



(48) 



This will be proven in the appendix. We stress that this spin 
contribution, apart from the spin quantum number, only de- 
pends on L — y . 

The quantity ( l44b can now be calculated in analogy to the 
case without spin^. Starting from equation ( l43T l, one employs 
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the expressions for AS from (HTI) and for w^™(s,u), the sum 
rule from Ref. 8 and the survival probability p (t), modified 



by replacing t with {t 
yields 



- I^^i {la - as in Ref. ^. This 



«2 m'^ 



{2s+l)!i 




E«(.,(rif*exp(-|)j£ 



, 1=1 
k 



AE 



2s +1 



N1+N2 l^i\Ni+N2j \2s+l 



L (-l)^A^(v) 

.L-V=k 



(49) 



The integrals over s and u were calculated in Ref. |^ and the 
sum over v can be carried out with the recursion formula^ 



V.L 



(50) 



2, 



where ^ — 1, if time reversal symmetry is present and p 
if time reversal symmetry is broken. 

Finally, using these results in the case of time-reversal in- 
variance, we obtain for the full transmission matrix, including 
also the diagonal part, 

2s+l (2s+l)2 



r-nd 



«2.«, A^i+A^2 {2s+l)iNi+N2)-l 
in the case of half-integer s, and 



2s+l 



Ni 



{2s+lf 



(51) 



(52) 



-N2 {2s+l){Ni+N2) + l 

if s is integer For s = 1/2 the result ( fSTT ) is identical with 
the one obtained using Random Matrix Theory, in the circular 
symplectic ensemble'. 

These findings can now be compared with the respective re- 
sults when time-reversal is absent, thus revealing the behavior 
of the transmission under a breaking of time-reversal by, e.g., 
turning on a magnetic field. In that case p = 2 so that the term 
(ISOl l vanishes, implying via ( |49] l that only the diagonal contri- 
bution survives. The difference AT — T 'P^'^ — T '^P"^' of the 
transmission coefficients therefore is 

NiN2{2s^\) 



AT w 

(M+A^2)((2.9+1)(M+A^2)- 

in the case of half-integer s, and 



1) 



AT 



-MA^2(2.?+1) 



1 



(53) 



(54) 



(M+A?2)((2^+1)(M+A^2)- 

if s is integer From these expressions one immediately con- 
cludes that the transmission (i.e., conductivity) is enhanced at 
zero magnetic field (when time reversal symmetry is restored), 
if the spin is half-integer; thus weak anti-localization occurs. 
To the contrary, integer spin would lead to weak localization. 
The latter had previously been obtained in semiclassical stud- 
ies where the spin had been neglected^. The only semiclassi- 
cal derivation of weak anti-localization so far'^, however, was 
restricted to the one-loop contribution and employed asymp- 
totics for large Ni,N2- 



r 



IV. SHOT NOISE 



The techniques developed above can be applied to a number 
of further problems arising in the context of ballistic transport 
through chaotic mesoscopic cavities. As a first example we 
consider shot noise. To this end one needs to compute the 
energy-averaged Fano-factor F, defined as 



F := 



(Tr(rr' -rr'rr')) 



AE 



(Tr(rr+)) 



(55) 



AE 



in terms of the transmission matrix T = S^^. The denominator 
has been dealt with above, and the spin-independent contribu- 
tion to 



(56) 



was calculated semiclassically in Ref. I?. We are hence left 
with the task of determining the spin contribution to ( l56t . Re- 
ferring to the semiclassical representation (VT\ one immedi- 
ately realizes that a four-fold sum over classical trajectories 
emerges. In addition to the case covered in Ref. JJ^ each term 
in this sum acquires an additional factor of 



Tr 



(d\DuDId„ 



(57) 



in which the indices label the trajectories involved. The diag- 
onal contribution to the four-fold sum occurs with s ~ u and 
V = w, or with s — w and u = v. In both cases unitarity implies 



Tr D]D„DlD 



:25+l 



(58) 



Beyond this one has to consider the encounter of four trajecto- 
ries. For the first time this has been done in quantum graphs'^, 
and has later been extended in Ref.[3- Following the method 
of these papers, every trajectory consists of two parts, labeled 
by 1 and 2. Approximately, one then has si — w\, mi = vi, 
52 = U2 and V2 = 1^2 . Thus 



Tr 



2s+l. (59) 
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Following further the calculation of the Fano factor in Ref.lli 
we obtain 



N1N2 



iNi+N2y 



(60) 



for Ni,N2 ^ 1. This result coincides with the respective out- 
come of a random matrix calculation in the symplectic ensem- 
ble 



,1.16 



V. CONDUCTANCE FLUCTUATIONS 

Universality of conductance fluctuations is often character- 
ized in terms of the energy-averaged variance of Tr(rr^). 
Instead of this quantity, the energy-averaged covariance of 
Tt{R"R"^), where n = 1,2 labels the leads, can also be con- 
sidered, see Ref.[ll|for details. Our calculations are based on 
the first paper of Ref. [ll], whose method can still be applied 
when the Ehrenfest time is much smaller than the dwell time; 
this condition is fulfilled in the semiclassical limit considered 
here. 

The calculation of the variances again involves four-fold 
sums over trajectories, in which the spin contribution occurs 
in terms of the factors 



Tr(D,Dl^Tr(D,Di 



(61) 



Switching off the spin-orbit interaction while preserving the 
presence of spin s, one obtains 

Tr Tr (d^) = (2. + 1 )2 . (62) 

In the presence of spin-orbit interaction one must examine the 
trajectories involved more closely. Here we again consider 
the case Ni,N2^ 1. The trajectories are divided into three 
parts labeled by 1, 2 and 3, and the relations si — ui, S2 — 

V2, S3 — M3, Vi = Wl, U2 = W2, V3 = W3 Or Si = Ml, S2 — V2, 

S3 = M3, vi = Wl, M2 = W2, V3 = W3 hold approximately. Here 
an over-bar indicates that these pieces are traversed in reverse 
direction. In the first case this yields 

Tr (DsDI) Tr (£),Z)+ ) « Tr (As^oJ,) Tr (z)„,d;) , (63) 

whereas in the second case 

Tr ) Tr (d„d1,) « Tr ) ' . (64) 

After an average over SU(2), very much alike in the main part 
of this work, we obtain for the first case22i2i 



75(7(2) JSU(2) 

and for the second case'- 



dgadgbTr (n, (gagt) ) Tr (n, (gbgl 



I I 

JSU(2) JSi 



SU{2) JSU{2) 



?aagb 



Tr 71,. ga 



1 , 

(65) 
(66) 



We follow Ref. [TT| further and finally observe that, with 
M ,A^2 ^ 1, the energy-averaged variance of Tr(rr^) reads 



when the spin-orbit interaction is switched off, and 

in the presence of spin-orbit interaction. Again, this finding 
is in accordance with the respective result in the symplectic 
ensemble ofRMTii^. 



VI. SUMMARY AND CONCLUSIONS 

We considered the semiclassical description of ballistic 
transport through chaotic mesoscopic cavities in the presence 
of spin-orbit interactions. Our focus was the calculation of 
transmission coefficients. Here the principal task was to ver- 
ify the effect of weak anti-localization in the form predicted 
by RMT. 

Working within the framework of the Landauer formalism, 
our starting point was a semiclassical representation of Green 
functions for Hamiltonians that contain a spin-orbit interac- 
tion. Transmission coefficients then require the evaluation of 
double sums over classical trajectories. The principal diffi- 
culty presented by such expressions is to get hold of the inter- 
ferences thus occurring. This can be overcome successfully 
by exploiting the Sieber-Richter method, originally developed 
to perform analogous calculations in the context of spectral 
fluctuations in classically chaotic quantum systems. 

We attacked the problem using the two established vari- 
ants of the Sieber-Richter method: the configuration-space ap- 
proach for the leading order, and the phase-space approach for 
the remaining contributions. In the first case a key input was a 
classical sum rule encoding an ergodic (and mixing) behavior 
of the combined classical spin-orbit dynamics. Essential to 
the success of the phase-space approach was a calculation of 
the spin contribution to pairs of classical trajectories that are 
grouped together pairwise according to the structure of their 
almost self-encounters. This led to the central result given in 
Eq. ( |48] l. The sign appearing points to the essential difference 
between the effects of half-integer spin as opposed to integer 
spin (including spin zero). This difference was then identi- 
fied as responsible for weak anti-localization or localization, 
respectively, to occur. We finaUy showed how our approach 
generalizes to semiclassical descriptions of shot noise and of 
universal conductance fluctuations. 



Appendix A: PROOF OF THE RELATION @8) 

We will show the validity of Eq. (|48] l by induction with 
respect to the number n of 2-encounters of two trajectories 
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"i^'i ■ The proof is based on the relations 

/ dgTr{%,{xgyg))=^-^Tx{%,{xy-')) (Al) 
ist/(2) 2^ + 1 ' 

and 

JSU{2) JSU{2) 

= ^2^^'^'''^^'''''^^' ^^^^ 

valid for all w, x,y, z E SU (2). For finite groups analogous 
identities have been shown in Ref. |2ll their proofs can be di- 
rectly carried over to the present case. 
We now proceed in three steps: 

1. First consider the case n — 0, where y' — J- This also 
I 



means f\j = I and kj = j. Here we obtain 

= ... dgL...dg2 

JSU{2) JSU{2) 

xTr(jt, [gL-gisl-gl)) 
= 2s+l. (A3) 

2. We assume the validity of ( |48] ) for two trajectories y = 
{Ii,a,b,l2,c,d,l3) andy = {I4,a,b,l5,c,d,l6) as shown 
in Figure [3] Here stands for stretches of the tra- 
jectories Y and "/ containing an unspecified number of 
2-encounters. By assumption, the actual number of 2- 
encounters, where y differs from "/ is n. We show now 
that the relation ( |48] | is stiU vaUd, when we replace y 
with the trajectory y" = (Z4, a, c, Z5 , b, d,l(,). Thus y" dif- 
fers from jinn' — n+ I 2-encounters. Then 



MyY' = / •••/ dgadgbdgcdgd...Ti(nJgi^gdgcgi2gbgagiiglglgcgiigbglgl 

J SU (2) J SU (2) ^ ^ 

= / .../ dgxdgydg,...TT(nJgi,gxgt2gyghglglgz8i58z8l8i 

JSU{2) JSU{2) \ \ ^ '6 

= V~rT / ■ ■ ■ / dgadgbdgcdgd . . . Tr [71,. (gi^gdgcghgbgagkg]^glglg\glglg]^ ) 

■^S -\- V J SU (2) J SU (2) ^ ^ ■' 

(-1)2^- 



In the second step we substituted gjgc ~ gx, gtgc = gz, 
and gijga = gy, and in the third one we used Eq. ( lAll i. 
In the fourth step we undid the substitution. This calcu- 
lation proves that changing the number of 2-encounters, 
in which y and y differ, by one indeed contributes a fac- 
tor of (-1)2-7(2.9 + 1). 

3. We assume the validity of the rela- 
tion ( |48] l for the two trajectories y = 

(li,a\,b\,l2,a2,b2,h,ci,d\MiC2,d2M) and 
y = {I6,ai,bi,h,a2,b2,k,ci,di,l9,c2,d2,lw) as 
shown in Figure |4] Again we assume that the 



number of 2-encounters, where y differs from 
is n. We then show that the relation ( l48T l is un- 
changed under a replacement of y with the trajectory 
y = {le,ai,di,lg,C2,b2,k,ci,bi,h,a2,d2,liQ). Notice 
that y cannot be constructed by applying the procedure 
of 2. twice: here the stretches /g, Ij and of y are tra- 
versed in parallel direction, whereas in 2. the stretches 
I4 and ^6 of y are traversed in anti-parallel direction. 
A calculation similar to (IA4l i. with the substitutions 
8dj8cj = 8xj, 8hj8lj = 8zj, 8bj8aj = 8yj U ^ {1,2}), 
then yields 
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Y/' 



SU{2) JSV{2) 



.Tr 



("i {gligchgciEkgch get gkgb2ga2gl2gbi gcngh gisWdi glggliSliSlgli gli gUligWhc 



SU{2) 

1 



SU{2) 



dgxi ■ • - Tr {ghgx2gkgxighgy2ghgy\ghg\gligzig]ggl2gl2gkgligligh^^^^ 



{2s+\Y jsu{2)"' jsu(2) 



dgx.dgy^dg^^dgy^ . . . Tr [n^ (^g\gl2g\gligXgl2s\aghgx2gkgxighgy2ghgyighg]f,gli 



1 



(A5) 



After these steps ( l48T l follows by induction because every tra- 
jectory y can be constructed successively out of y by using 
the procedures of 2. and 3. Every /-encounter that does not 
decompose into several encounters of a lower number of tra- 
jectories (see Figure 4 in Ref. 6 for an example) can be con- 
structed from 2-encounters in / — 1 steps. Every such step 



then brings out a factor of (— 1 )^*/ (2,? + 1 ) in My y , when this 
is constructed from My ^ = 2s + 1 . Thus, V encounters with al- 
together L stretches contribute a factor ((— 1)^ V (2s + 1))^ ^ , 
which completes the proof of ( |48t . 



* Electronic address: jens.bolte@uni-ulm.de 

^ Electronic address: DanieI.WaItner@physik.uni-regensburg.de' 

^ C.W.J. Beenakker, Random-matrix theory of quantum transport. 
Rev. Mod. Phys. 69 731 (1997) 

^ R. Landauer, Spatial variation of currents and fields due to lo- 
calized scatterers in metallic conduction, IBM J. Res. Develop. 1 
223 (1957); IBM J. Res. Develop. 32 306 (1988) 



H. U. Baranger, R. A. Jalabert, A. D. Stone, Quantum-chaotic 
scattering effects in semiconductor microstructures , Chaos 3 665 
(1993); Weak Localization and Integrability in Ballistic Cavities, 
Phys. Rev. Lett. 70 3876 (1993) 

M. Sieber, K. Richter, Correlations between Periodic Orbits and 
their Role in Spectral Statistics, Physica Scripta T90 128 (2001) 
^ S. Mtiller, S. Heusler, P. Braun, F. Haake, A. Altland, Semiclas- 



12 



sical Foundation of Universality in Quantum Chaos, Phys. Rev. 
Lett. 93 014103 (2004) 

* S. Heusler, S. Muller, P. Braun, F. Haake, A. Altland, Periodic- 
orbit theory of universality in quantum chaos, Phys. Rev. E 72 
046207 (2005) 

^ S. Heusler, S. Miiller, A. Altland, P. Braun, F. Haake, Periodic- 
Orbit Theory of Level Correlations, Phys. Rev. Lett. 98 044103 
(2007) 

^ K. Richter, M. Sieber, Semiclassical Theory of Chaotic Quantum 

Transport, Phys. Rev. Lett. 89 206801 (2002) 
^ S. Heusler, S. Muller, P. Braun, F. Haake, Semiclassical Theory of 

Chaotic Conductors, Phys. Rev. Lett. 96 066804 (2006) 
" S. Muller, S. Heusler, P. Braun, F. Haake, Semiclassical approach 

to chaotic quantum transport. New J. Phys. 9 12 (2007) 
' PW. Brouwer, S. Rahav, A semiclassical theory of the Ehrenfest- 

time dependence of quantum transport in ballistic quantum dots, 

preprint, arXiv:cond-mat/05 12095 i'2 (2005); Phys. Rev. B 74 

075322 (2006) 

^ G. Berkolaiko, J.M. Harrison, M. Novaes, Full counting statis- 
tics of chaotic cavities from classical action correlations, preprint, 
^Xiv:cond-mat/0703803 ^ (2007) 

^ H. Schanz, M. Puhlmann, T. Geisel, Shot Noise in Chaotic Cavi- 
ties from Action Correlations, Phys. Rev. Lett. 91 134101 (2003) 

* S. Heusler, S. Muller, P. Braun, F. Haake, Semiclassical Predic- 
tion for Shot Noise in Chaotic Cavities, J. Phys. A: Math. Gen. 39 
L159 (2006) 

^ I. Zutic, J. Fabian, S. Das Sarma, Spintronics: Fundamentals and 
applications. Rev. Mod. Phys. 76 323 (2004) 



D. V. Savin, H.-J. Sommers, Shot noise in chaotic cavities with an 
arbitrary number of open channels, Phys. Rev. B 73 081 307 (R) 
(2006) 

O. Zaitsev, D. Frustaglia, K. Richter, Role of orbital dynamics in 
spin relaxation and weak antilocalization in quantum dots, Phys. 
Rev. Lett. 94 026809 (2005); Semiclassical theory of weak an- 
tilocalization and spin relaxation in ballistic quantum dots, Phys. 
Rev. B 72 155325 (2005) 

J. Bolte, S. Keppeler, A semiclassical approach to the Dirac equa- 
tion, Ann. Phys. (NY) 274 125 (1999) 

J. Bolte, S. Keppeler, Semiclassical form factor for chaotic sys- 
tems with spin 1/2, J. Phys. A: Math. Gen. 32 8863 (1999) 
J. Bolte, J. Harrison, The spin contribution to the form factor of 
quantum graphs, J. Phys. A: Math. Gen. 36 L433 (2003) 
J. Bolte, J. Harrison, The spectral form factor for quantum graphs 
with spin-orbit coupling, in: G. Berkolaiko, R. Carlson, S.A. 
Fulling, and P. Kuchment (eds.): Quantum Graphs and Their 
Applications, Contemporary Mathematics, Volume 415, pp. 51 
(AMS 2006) 

D.S. Fisher, PA. Lee, Relation between conductivity and trans- 
mission matrix, Phys. Rev. B 23 6851 (1981) 
H.U. Baranger, A.D. Stone, Electrical linear-response theory in 
an arbitrary magnetic field: A new Fermi-surface formation, 
Phys. Rev. B, 40 8169 (1989) 

J.U. Nockel, A.D. Stone, H.U. Baranger, Adiabatic turn-on and 
the asymptotic limit in linear-response theory for open systems, 
Phys. Rev. B 48 17569 (1993) 



